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SUMMARY 


This paper investigates the stability of the negative feedback interconnection 
of two positive-real systems which have poles in the closed left-half of the complex 
plane. A new definition of marginally strict positive real systems is introduced, 
and sufficient conditions are obtained for the stability of the feedback 
interconnection of such systems, using the Lyapunov method. The conditions obtained 
have direct applications to dynamic dissipative controllers for flexible spacecraft, 
and are the least restrictive ones published to date. 

I. INTRODUCTION 

A well known result in passivity theory [1] is the Passivity Theorem, which 
states that a stable passive dynamic system can be stabilized by a stable strictly 
passive controller. For linear, time-invariant (LTT) systems, passivity is 
equivalent to positive realness. However, the relationship between strict passivity 
and strict positive realness is somewhat more complicated because there are several 
definitions of strict positive realness. Strict passivity is equivalent to the 
strongest definition of strict positive realness (SPR) [2] for LTI systems. (We 
shall refer to this definition simply as "strict passivity", even for the LTI case). 

The Passivity Theorem, which is for general (possibly nonlinear) systems, can be 
directly applied to LTI systems; i.e., the feedback interconnection of a stable PR 
system and a strictly passive system is stable. However, the requirement of strict 
passivity is too stringent, as it includes only systems with a relative degree of 
zero. Another definition of strict positive realness was introduced in [2] for 
scalar systems, and was further investigated in [3, 4, 5] for multivariable systems. 

This definition (referred to as "strong SPR") is weaker than strict passivity. 

In [6], an even less stringent definition of SPR, termed "weak SPR", was 
investigated, and it was stated that the feedback interconnection of a stable PR 
system and a weak SPR system is stable. That is, the strict passivity requirement 
was replaced by the "weak SPR" requirement. This was a significant improvement. 



since weak SPR systems can include strictly proper systems, while strict passivity 

requires the relative degree to be zero. However, in all the studies in the 

published literature, both weak and strong definitions of SPR require the systems to 

be stable (i.e., all poles are required to be in the open left-half of the complex 

plane). In this paper, we shall remove this restriction. 

Let G(s) denote an mxm matrix whose elements are proper rational functions of 

the complex variable "s". Let the conjugate-transpose of a complex matrix T be 
* 

denoted by T . 

Definition 1: An mxm rational matrix G(s) is said to be positive real (PR) if 

(i) all elements of G(s) are analytic in Re[s] > 0; 

3 (£ 

(ii) G(s) + G (s) > 0 in Re[s] > 0 , or equivalently, 

(iia) poles on the imaginary axis are simple and have Hermitian, 

nonnegative- definite residues, and 

(iib) G(jto) + G (jto) ^ 0 for co € (-°°,°°) 

Given below are some definitions of strictly positive real systems. Definition 
2, which represents the specialization to LTT systems of the general definition of 
strict passivity, is the strongest definition of strict positive realness. 

Definition 2: An mxm rational matrix G(s) is said to be strictly passive if 

(i) all elements of G(s) are analytic in Re[s] > 0; 

(ii) there exists an e > 0 such that 

ifj 

G(jto) + G (jco) > el for to e (-°°,°°) 

Definition 3: An mxm rational matrix G(s) is said to be strictly positive real in the 
strong sense (strong SPR, or SSPR) if G(s-e) is PR for some e > 0; that is, if 

(i) all elements of G(s) are analytic in Re[s] > 0; 

(ii) G(jto) + G (jco) > 0 for to e (- 00 , 00 ) 

(iii) Z= G(°o)+G T (oo) > 0 

(iv) ^[GCjto) + G (jto)] > 0 if Zis singular 
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Definition 4: An mxm rational matrix G(s) is said to be strictly positive real in the 
weak sense (weak SPR, or WSPR) if 

(i) all elements of G(s) are analytic in Re[s] > 0; 

(ii) G(j(o) + G (jco) > 0 for to e (-«,«>) 

Note that Definition 2 requires that Z= G(°°)+G T (°°) to be positive definite; 
i.e., the system must have a relative degree of zero. This requirement makes the 
definition of strictly passive systems too restrictive. Definition 3 (SSPR) can 
include certain strictly proper systems which satisfy additional conditions (iii) and 
(iv). Definition 4 (WSPR) does not require these additional conditions, and is 
therefore less restrictive than Definition 3. However, all the definitions (2-4) of 
SPR require the system to be stable. 

In this paper, we go one step further, and allow the system to have poles on the 
imaginary axis. The significance of this is that many physical systems can now be 
included in this (much larger) class of systems. 

Definition 5: An mxm rational matrix G(s) is said to be marginally strictly positive 
real (MSPR) if it is positive real, and 

jjc 

G(jco) + G (jco) > 0 for co € (-°°,«0 

Definition 5 of MSPR differs from Definition 1 (PR) because the frequency domain 

inequality (>) has been replaced by the strict inequality (>). The difference 

between Definitions 4 and 5 is that Definition 5 allows G(s) to have poles on the 

imaginary axis. This is an important difference because many real-life systems 

contain pure integrators and oscillators, which are permitted under Definition 5, but 

y 8 s 

not under Definitions 2, 3, and 4. For example, let G(s) = -h- — — +H(s) , where y 

s s 2 +0) 2 
o 

and 6 are real non-negative scalars and H(s) is weak SPR. Then G(s) is marginally 
SPR. 
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Suppose [A,B,C,D] is a minimal realization of a rational matrix M(s). M(s) (or 

[ A,B,C,D ]) is said to be minimum-phase if its transmission zeros are confined to the 

[ s / -A B1 

q D I can drop below its normal value only 

for values of s in the OLHP. 


n. PROPERTIES OF MSPR SYSTEMS 

Suppose G(s) is positive real and has all poles (i.e., the eigenvalues of the 
system matrix of its minimal realization) in the closed left-half plane (CLHP). 
Following [7], G(s) can be written as: 

G(s) = G^s) + G 2 (s) (1) 


where G^s) has purely imaginary poles, and G 2 (s) has poles only in the open 
left-half plane (OLHP). Furthermore, G^s) is of the form: 

p 


G,(s) 


a s +B 

_» + y_j — ^ 

S ,^-S 2 + O ) 2 


i = 1 


( 2 ) 


where a., and (3. are mxm real matrices, and G).>0, i=l,2,...,p (to.# (0. for i*j). 

Some remarks regarding the nature of the poles on the imaginary axis are in 
order. The poles and zeros considered here are in the Smith-McMillan sense [8]; 
i.e., there can be more than one pole at a given location, without it being 
considered a "repeated" pole. In particular, using standard results in matrix 
fraction descriptions [8], it can be shown that the McMillan degree (i.e., the 
minimal order of a state space representation) of the term: [a Q /s] is equal to p(a Q ), 
where p(.) denotes the rank. That is, there are v=p(a Q ) simple poles at s=0. 

Suppose the McMillan degree of the term: [as+B ]/(s + 0 ) ) is 2k , where k <m. 

i i i i i 

Then this term has k. simple poles (each) at s=jco. and s=-j(t).. 

The following results state that G(s) is PR (respectively, MSPR) iff its stable 
part is PR (resp. WSPR). 
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Lemma 1. G(s) is PR (resp. MSPR) iff all of the following hold; 

(i) G 2 (s) is PR (resp. WSPR); 

(ii) a.= £ 0, i =0, 1, 2, ...p; 

(iii) P,= i =1 ,2, ...p. 

Proof. From the requirement that the residues at the imaginary- axis poles be 
nonnegative-definite, we get (ii) and (iii) (See [9]). Therefore, we have: 

G^co) + G*(jo)) = 0 , and G(jco) + G*(j(fl) = G 2 (jco) + G*(jco), which is positive 
semi-definite (resp. positive definite) for all real to iff G 2 (s) is PR (resp. WSPR). ■ 

* 

If G(s) is MSPR, the degree of [G(s)+G (s)] is generally less than the degree 
of G(s); also, G 2 (s) is stable, and both G(s) and G 2 (s) are minimum-phase. 

We shall next consider a minimal realization of MSPR transfer functions. 


III. MINIMAL REALIZATION OF MSPR SYSTEMS 

Consider the realization of Gjfs). For the term [(Xg/s], since a Q is symmetric 
and non-negative definite (Lemma 1), there exists an mxm real orthogonal 
transformation matrix T which diagonalizes it, i.e., 

T t cc o T = diag[X I ,X 2 ,...^ v ,0,...,0] (3) 

where X. are positive scalars. Let A denote diag[A, ,X .....Xy]. A minimal 
realization (of order v) of T T [a o /s]T is: 0], where 


•w 


®o = Pv °* 



( 4 ) 


where 0 V and I denote vxv null and identity matrices. Therefore, a minimal 
(v'b-order) realization of [cc q /s] is given by [„/£, 0, where 
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( 5 ) 


GB= 0 T ; 
o o 


<?=T<? 
o o 


Considering the term: [a.s + p . ]/(s 2 +co 2 ), if its McMillan degree is 2k. (k.<m), 


a minimal realization is given by [7]: [A,0,0 r ,O] where A e R 


2kX2k 

> i. 


where 


A = diag(A ] ,fl 2 ,A k ) 


V R '»] 


( 6 ) 


(7) 


2k Xm 

and 0e R i . Then a minimal realization of G^s) is given by: [A^Bj.C .0], where 
A= diag(^, A, A 2 ,....,A); B ] = [0 q t , 0^...,0jf-, C = [$0 T ,...,0 T ]; Ae 


R n i“i, B,e IR 11 ™, Ce R mxn i (n = v+2? k.) 

1 = 1 * 


1 


1 


1 


Let [A 2 ,B 2 ,C 2 ,D] be a minimal realization G 2 (s), where A 2 e IR n 2 xn 2 , B 2 e (R n 2 xm , 

C 2 e !R mxn 2 , and D e IR mxm (n 2 =n - n^. Let A = diag(A } , A 2 ); B = [B^,B 2 ] T ; C = 

[C,CJ. Then A e (R nxn , B e IR nxm C e R mxn , and D e IR mxm . Since G (s) and G (s) 
12 12 

have no poles in common, [A,B,C,D] is a minimal realization of G(s). 


IV. CHARACTERIZATION OF MSPR SYSTEMS 

We next present the state-space characterization of MSPR systems, which is an 
extension of the Kalman- Yacubovich lemma for WSPR systems, that was proved in [6]. 


Lemma 2. If G(s) is MSPR, there exist real matrices: P = P 1 > 0, P e R nxn , _^?g R kxn 2 , 


W g R^" 1 (where k > m), such that 


A t P + PA = -L t L 
C = B t P +W t L 
W T W = D + D t 


6 


( 8 ) 

(9) 

( 10 ) 



( 11 ) 


L = »>,. -<EJ 

1 2 

where [A 2 ,B 2 ,-^fW] is minimal and minimum-phase. 

Proof. G(s) has the form given by Eqs. (1) and (2). Since G(s) is MSPR, conditions 
(i)-(iii) of Lemma 1 hold. Considering G^s) in Eq. (2), it consists of (p+1) 
transfer functions in parallel. A minimal realization (Eqs. 4 and 5) of the first 
term, [a Q /s], is: <^,0]. Letting fI o = A, it can be verified that the equations 

are satisfied: 


*As n + n yc — o 
0 0 0 0 


( 12 ) 


n n 0= 

0 0 0 


(13) 


A minimal realization (Eqs. 6, 7) of the ith component of the second term, 
[a.s+p.]/(s 2 +(0 2 ), is given by ,6fl T ,0], Letting II. = I . , it can be seen that 

i i 1 ill i 2K 

i 

this realization satisfies: 


yc n + n o 

i i i i 

(14) 

n SB= SB 

i i i 

(15) 

Finally, following [6,7], since G 2 (s) is WSPR, there exist P 2 = P 2 > 0, P 2 e 
e tR kxm , Jfe. R 1 ™^ (where k > m) such that 

R n 2 Xn 2 , W 

A P + P A = -Jf ¥ 
2 2 2 2 

(16) 

c = B^P +W T ^f 
2 2 2 

(17) 

W T W = D + D t 

(18) 


where [A ,B ..^W] is minimal and minimum-phase. 
Defining 


P = diag(n o ,n i ,...,IT p ) (19) 
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and 


P = diag(P if P 2 ) (20) 

we have the required result. ■ 

V. STABILITY OF FEEDBACK INTERCONNECTION 

Consider the system in Figure 1, where G(s) and H(s) are mxm proper rational 
matrices. This system is said to be stable if its state-space realization consisting 
of individual minimal realizations of G(s) and H(s), is asymptotically stable. We 
have the following stability result. 

Theorem 1. The negative feedback interconnection of G(s) and H(s) is stable if all of 
the following conditions are satisfied: 

(i) G(s) is MS PR; 

(ii) H(s) is PR; 

(iii) None of the joa-axis poles of G(s) is a transmission zero of H(s). 

Proof. Let [A,B,C,D] denote the minimal realization of G(s) described in Section III, 
and let x be the corresponding state vector of order n. Let n 2 denote the number of 
poles of the stable part G 2 (s) of G(s), and let [A 2 ,B 2 ,C 2 ,D] denote its minimal 

A A A A A 

realization. Let [A,B,C,D] denote a minimal realization of H(s), and let x denote 
the corresponding n order state vector. Since G(s) is MSPR, from Lemma 2, there 
exist matrices P = P T > 0, P e R nxn , R mxn 2, W e IR mxm , such that Eqs. (8)-( 1 1) are 
satisfied, and [A ,B 2 »-^fW] is minimal and minimum-phase. 

AA A AA A AA A A 

Since H(s) is PR, there exist matrices P = P T > 0, P e R nxn , L e R*" 1 , W e R 10 ^ 1 

such that [7] 


A x A A A Ax A 

a t p + PA = -l t l 

(21) 

A AtA A tA 

C = B P +W L 

(22) 

A t' A A At 1 

W‘W = D + D 
8 

(23) 


Consider the candidate Lyapunov function: 

A x Ax A A 

V(x,x) = x T Px + x T Px 

Proceeding as in the proof of Theorem 1 in [6] we have: 

. T T AxA AxA 

V = 2u y - z z + 2u y -z z 

where 

z = Lx + Wu = fx + Wu 
2 

A AA A A 

z = Lx + Wu 

A A 

Since u = y and y = -u, 

V=-zz-zz<-zz 

i.e., V is negative semi-definite. V=0 implies z=0. However, z(t) is the output 
the system: F(s) = W+L(sI-A) 'B = W+^sI-A )’ ! B , which is minimum-phase, 
mput u(t) that results in z(t) = 0 must have the form: u(t) =£ u.t j e b j , where 
C.’s are zeros of F(s), and Re[^ ] < 0. Therefore, u(t) -» 0 exponentially; i.e., 

A 

y(t)-»0 exponentially. As a result, y(t) will consist of i) exponentially decaying 
terms corresponding to the zeros of F(s) as well as stable poles of G(s), and ii) 
persistent terms such as corresponding to unstable poles at s=jco (including 

at s=0) of G(s). Since jto. are not the transmission zeros of H(s), this would imply 
that y(t) [=- u(t)] will contain persistent terms such as e^i 1 . However, this 
contradicts the fact proved previously that u(t) decays exponentially. Therefore, 
y(t) can consist only of exponentially decaying terms, i.e., y(t)->0 exponentially. 

A A A A A 

Because of the minimality of [A,B,C,D] and [A,B,C,D], this implies x(t)->0, x(t)-*0 
exponentially. Using LaSalle’s invariance theorem [10], the system as asymptotically 
stable. ■ 

It should be noted that, in Theorem 1, G(s) and H(s) are completely 
interchangeable. 

The following corollaries are an immediate consequence of Theorem 1. 
Corollary 1.1 The negative feedback interconnection of G(s) and H(s) is stable if 
G(s) is WSPR and H(s) is PR. 


(24) 

(25) 

(26) 

(27) 

(28) 
of 

Every 
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Corollary 1.1 is the same stability result which was given in [6]. 

Corollary 1.2 The negative feedback interconnection of G(s) and H(s) is stable if 
both G(s) and H(s) are MSPR. 

Example: Consider the rotational motion of a flexible spacecraft with m torque 
actuators and m collocated attitude sensors (m>3). Assume that there is at least one 
torque actuator for each (orthogonal) axis of rotation. The transfer function from 
the torque input to the attitude (position) output, y , is given by: 

p 

G(s) = G'(s)/s (29) 

where 

p 

OC i * cx s 

G'(s) = — + ) -r r (30) 

s *- — 's +2p co s+co 

i = l i 

where G'(s) is the transfer function from the torque input to the attitude rate 
y r (=y ); a,= oc^ > 0 (i =0, 1, 2, ...p), and a Q is a rank-3 matrix; (O.(>0) represent 
the natural frequencies, and p.>0 represents the inherent damping ratio, for the i 
elastic mode (i=l,2,..,p). It can be easily verified that G'(s) is PR, and 
therefore, from Theorem 1, it can be stabilized by any MSPR controller. Let <?(s) 
denote an mxm stable transfer function which has no transmission zeros on the 
imaginary axis, and suppose H(s) = [t?(s)/s] is MSPR. Then H(s) stabilizes G'(s). 
Examining the block diagram in Figure 2, <?(s) stabilizes G(s). In other words, a 
flexible spacecraft, which has zero-frequency rigid-body modes as well as damped or 
undamped elastic modes, is stabilized by the controller £?(s) which has the above 
properties. The stability does not depend on the number of elastic modes, or the 
parameter values, and is therefore robust. 

VI. CONCLUDING REMARKS 

The concept of marginally strictly positive real (MSPR) systems was introduced, 
which allows poles on the imaginary axis, and is therefore less restrictive than the 
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previous definitions of strict positive realness. A state-space characterization of 
MSPR systems was obtained, and it was proved that the negative feedback 
interconnection of an MSPR system and a positive real (PR) system, is asymptotically 
stable. The result significantly extends the previous passivity-based stability 
results for linear time-invariant systems. 
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Figure 1 . Negative Feedback Loop. 
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Figure 2. Flexible Structure Example. 
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